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Abstract —This paper presents a cooperative framework for fireworks algorithm (CoFFWA). A detailed analysis of existing fireworks 
algorithm (FWA) and its recently developed variants has revealed that (i) the selection strategy lead to the contribution of the firework 
with the best fitness (core firework) for the optimization overwhelms the contributions of the rest of fireworks (non-core fireworks) in the 
explosion operator, (a) the Gaussian mutation operator is not as effective as it is designed to be. To overcome these limitations, the 
CoFFWA is proposed, which can greatly enhance the exploitation ability of non-core fireworks by using independent selection operator 
and increase the exploration capacity by crowdness-avoiding cooperative strategy among the fireworks. Experimental results on the 
CEC2013 benchmark functions suggest that CoFFWA outperforms the state-of-the-art FWA variants, artificial bee colony, differential 
evolution, the standard particle swarm optimization (SPSO) in 2007 and the most recent SPSO in 2011 in term of convergence 
performance. 

Index Terms —Swarm Intelligence, Fireworks Algorithm, Cooperative Strategy, Explosion Amplitude 

- ♦ - 


1 Introduction 

In the past two decades, many stochastic and population-based 
swarm intelligence (SI) algorithms have been proposed. These 
algorithms have shown great success in dealing with optimization 
problems in various application fields. SI refers to the ability 
presented by the swarm behavior to solve problems through 
the direct or indirect interaction among the agents with their 
environment (TJ. Usually, SI system consists of a population of 
decentralized and simple natural or artificial agents, the iteration 
among the agents lead to the emergence of intelligent ability j2j. 
By the observing and modeling of the cooperative swarm behavior 
of living creatures, artificial systems and the physical properties 
of non-living objects, researchers are trying to understand those 
mechanisms and use them to design new algorithms. 

The famous SI algorithms comprise ant colony optimization 
(ACO) by modeling the cooperative behavior among ant colony 
through environment to search for the shortest path from the 
colony to the food |[3j, and particle swarm optimization (PSO) by 
mimicking the aggregating motion of a flock of birds for searching 
for food Q. 

Recently, the great successes of ACO and PSO for solving 
optimization problems greatly push the SI developments forward. 
Algorithms inspired by collective biologic behavior of bee, glow¬ 
worm, fish school, firefly, cuckoo, krill herd @ 0 0 H) @ m 
© ITS) p~3) | l4j fl5) p~6) and collective artificial system are 
proposed [ 17 |18) |19| . See |S0) for a detailed review of recently 
developed computational intelligence algorithms. 
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(a) fireworks explosion (b) solution search 


Fig. 1. Comparison between fireworks explosion and solution search for 
optimization problems. 

Inspired by the phenomenon of fireworks explosion in 
the night sky, a new SI algorithm called fireworks algorithm 
(FWA) has been proposed recently. In FWA, each firework 
in the fireworks population performs a local search around its 
position, and cooperates with each other for the global search. 
The explosion of fireworks can be seen as the search of solutions 
in the feasible range (see Fig. |T|). Since FWA has shown success 
in solving practical optimization problems, it is believed that FWA 
is efficient and worth for further study. 

1.1 Related work 

Related work based on FWA can be grouped into three categories, 
theory analysis, algorithm developments and applications. 

1.1.1 Theory analysis 

In (2l| , Liu presented a theory analysis for FWA, and proved 
that FWA is an absorbing Markov stochastic process. Moreover, 
the convergence property and time complexity of FWA was also 
discussed. 
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1.1.2 Algorithm developments 

Most studies focussed on algorithmic improvements of single 
objective FWA developments, multi-objective FWA developments 
and parallel FWA implementations. 

In |22| |, Pei et ai. investigated the influences of approximation 
landscape approaches on accelerating FWA by introducing the 
elite point of the approximation landscape into the fireworks 
swarm. Additionally, Liu et ai. (23) investigated the influences 
on FWA when using new methods for computing the explosion 
amplitude and the number of explosion sparks while still main¬ 
taining the core idea that fireworks with smaller fitness will have 
more sparks and smaller explosion amplitudes. 

The first comprehensive study of the operators of FWA can be 
found in |24| , where Zheng S. et al. proposed the enhanced fire¬ 
works algorithm (EFWA), which incorporates five modifications 
compared to conventional FWA: ( i ) a new minimal explosion am¬ 
plitude check, ( ii ) a new operator for generating explosion sparks, 
(Hi) a new mapping strategy for sparks which are out of the search 
space, (iv) a new operator for generating Gaussian sparks, and ( v) 
a new operator for selecting the population for the next iteration. 
The limitations in conventional FWA are presented and several 
improvements are proposed to deal with these limitations. 

Based on the work of EFWA, Zheng S. et at. proposed the 
dynamic search fireworks algorithm (dynFWA) ]25) . In dynFWA, 
the firework with minimal fitness in each iteration is called core 
firework and uses the dynamic explosion amplitude strategy while 
for the rest of fireworks, they use the same strategies as in EFWA 
(see |3.2.l] for the details of dynFWA). In addition, Li ]26| proposed 
an adaptive version of FWA (AFWA), in which the explosion 
amplitude of the firework with minimal fitness is calculated as the 
infinite norm distance between the firework and a certain candidate 
from the explosion sparks (see |3.2.2| for the details of AFWA). 

Additionally, a number of hybrid algorithms between FWA 
and other algorithms have been proposed. Gao et al. (27) proposed 
cultural fireworks algorithm which is the hybrid between cultural 
algorithm and FWA. Zheng Y. et al. (28| proposed a hybrid algo¬ 
rithm FWA-DE of FWA and differential evolution (DE). Zhang et 
al. (29) proposed a hybrid biogeography-based optimization with 
fireworks algorithm while Yu et al. (30) also presented a hybrid 
algorithm between DE with FWA. 

For multi-objective FWA (MOFWA), Zheng Y.et al. proposed 
the framework of MOFWA which uses the fitness assignment 
strategy from |3l( and minimal pairwise distance metrics in (32) . 

For Parallel FWA implementation, Ding et al. proposed GPU- 
FWA (33j|, a GPU implementation which maintains a high speedup 
value. 

1.1.3 Applications 

For the practical applications, FWA has been applied in FIR and 
HR digital filters design (27) , the calculation of non-negative 
matrix factorization (NMF) (34)-|36) , spam detection 1371, image 
recognition (38) , power system reconfiguration scheme [39|, mass 
minimisation of trusses (40] 6D- non-linear equation set (42) and 
0/1 knapsack problems |43J]. 

All above applications indicate that FWA is one promising SI 
algorithm and has a profound value and wide future. 

1.2 Contributions 

The main contributions of this paper are: ( i ) an evaluation of 
the cooperative strategies among the fireworks of EFWA, dyn¬ 


FWA and AFWA is conducted, (ii) the proposal of the coop¬ 
erative framework for FWA with independent selection method 
and crowdness-avoiding cooperative strategy. In the cooperative 
framework, the independent selection method will ensure the 
information inheritance and improve the local search ability for 
each firework while the cooperative strategy can enhance the 
global search ability of the fireworks swarm. (Hi) a comparison of 
the two different explosion amplitude strategies in dynFWA and 
AFWA, as well as the relationship between them are presented. 

The remainder of this paper is organized as follows: Sec¬ 
tion [2] introduces the framework of FWA and EFWA, and gives 
a comparison between FWA with other meta-heuristic algorithms. 
Section [3] describes two state-of-the-art improvement works on 
FWA, i.e dynFWA and AFWA, which focus on the explosion 
amplitude strategies. In Section [4| we present a comprehensive 
analysis of cooperative strategies in conventional FWA framework. 
To overcome the limitations, the new cooperative framework of 
FWA is finally proposed. To validate the ideas, several experiments 
are conducted in Section [5] and [6] Finally, concluding remarks are 
given in Section [7] 


2 The Framework of FWA and EFWA 
2.1 General Framework 

A principal FWA works as follows: At first, N fireworks are 
initialized randomly, and their quality (i.e., fitness) is evaluated 
in order to determine their explosion amplitude (i.e., the explosion 
range) and the number of sparks for each firework. Subsequently, 
the fireworks explode and generate different types of sparks 
within their local space. In order to ensure diversity and balance 
the global and local search, the explosion amplitude and the 
population of the newly generated explosion sparks differ among 
the fireworks. A firework with better fitness can generate a larger 
population of explosion sparks within a smaller range , i.e., with a 
small explosion amplitude. Contrary, fireworks with lower fitness 
can only generate a smaller population within a larger range, i.e., 
with higher explosion amplitude. This technique allows to balance 
between exploration and exploitation capabilities of the algorithm. 
Exploration refers to the ability of the algorithm to explore various 
regions of the search space in order to locate promising good solu¬ 
tions, while exploitation refers to the ability to conduct a thorough 
search within a smaller area recognized as promising in order to 
find the optimal solution (44) . Exploration is achieved by fireworks 
with large explosion amplitudes (i.e., lower fitness), since they 
have the capability to escape from local minima. Exploitation is 
achieved by fireworks with small explosion amplitudes (i.e, high 
fitness), since they reinforce the local search ability in promising 
areas. After the explosion, another type of sparks are generated 
based on a Gaussian mutation of randomly selected fireworks. 
The idea behind this is to further ensure diversity of the swarm. 
In order to improve readability we assign new notations to the 
two distinct types of sparks: “explosion sparks” are generated 
by the explosion operator, and “Gaussian sparks” are generated 
by Gaussian mutation operator. To retain the information of the 
fireworks swarm, and pass it to the next iteration, a subset of 
the whole population is selected at the end of each iteration. 
The algorithm continues until the terminal criterion is met. The 
framework of (E)FWA is given in Alg. [7] 
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Algorithm 1 - The general structure of (E)FWA. 

1: Initialize N fireworks and set the constant parameters 

2 : repeat 

3: // Explosion operator 

4: Calculate the explosion amplitudes 

5: Calculate the numbers of explosion sparks 

6: Generate “explosion sparks” (check if out-of-bounds) 

7: Evaluate the fitness 

8: // Gaussian mutation operator 

9: Generate “Gaussian sparks” (check if out-of-bounds) 

10: Evaluate the fitness 

11: // Selection strategy 

12: Select fireworks for the next iteration 

13: until termination criterion is met. 


2.2 Explosion Operator 

To perform the search for a firework, the points around the 
firework’s position are sampled with uniform distribution ( taking 
the inspiration from fireworks explosion process in the nature, see 
Fig. |T|). For this sampling operation, the sampling points number 
and the sampling range are two important factors which will 
influence the optimization results. The sampling points number 
is denoted as explosion sparks number and the sampling range is 
denoted explosion amplitude. 


2.2.1 Explosion amplitude and explosion sparks number 
Assume that the optimization problem / Q is a minimization 
problem, the fireworks number is N, the explosion amplitude A 
(Eq. Q) and the number of explosion sparks s (Eq. [2]) for each 
firework X, are calculated as follows. 


A i = A- 
Si = M e ■ 


f ( X i) Vmin A 0 
EZl (/ ( X i) - 2/min) + £ 
2/max H" £ 

E/i=l (^/max — f( X i )) + £ 


(1) 

( 2 ) 


where, y max = max(/(Xj)), 

2/min — min(/(Xj)), and A 
and M e are two constants to control the explosion amplitude 
and the number of explosion sparks, respectively, and e is the 
machine epsilon. Additionally, the number of sparks s, that can be 
generated by each firework is limited by the lower/upper bounds. 

Eq.[l]reveals that the explosion amplitude of the firework at the 
best location Xi, is usually very small [close to 0], the explosion 
sparks for Xi, will be located at the same location as Xb- As a 
result, these sparks cannot improve the location of X k . To over¬ 
come this problem, EFWA used the minimal explosion amplitude 
check strategy (MEACS) to bound the explosion amplitude A,; of 
firework Xi as follows. 


Ai = 


A ■ 

-^min 

A, 


if Ai A m j n , 
otherwise, 


(3) 


where, Aj" nin decreases non-linearly with increasing number of 
function evaluations such that 


Amin = A init - ^ init ^ final (4) 

A ma x 

1. In this paper, without loss of generality, the optimization problem / is 
assumed to be a minimization problem. 


where, t refers to the number of function evaluation at the 
beginning of the current iteration, and F max is the maximum 
number of evaluations. Aj n j t and Afj na i are the initial and final 
minimum explosion amplitude, respectively. 

2.2.2 Generating explosion sparks 

Now, each firework explodes and creates a different number of 
explosion sparks within a given range of its current location. For 
each of the explosion sparks of each firework X,. Algorithm [2] 
is performed once. 


Algorithm 2 - Generating one “explosion spark” in EFWA 
1: Initialize location of “explosion sparks”: Xi = Xi 
2: Set z k = round(rand(0,1)), k = 1,2,..., d. 

3: for each dimension of where z k == 1 do 
4: Calculate offset displacement: AX k = A,; x rand(—1,1) 

5: Xik — Hik A AX 

6: if Xik is out of bounds then 

7: randomly map Xik to the potential space 

8 : end if 

9: end for 


2.2.3 Mapping out-of-bounds sparks 

When the location of a new explosion spark exceeds the search 
range in dimension k, this spark will be mapped to another loca¬ 
tion in the search space with uniform distribution (in dimension k) 
according to x lk = X^ in + rand(0,1) * (X^ ax - X^J, here 
rand(a, b) denotes a random number from uniform distribution in 
[a,b\. 


2.3 Gaussian Mutation Operator 

In addition to the explosion sparks for local search, FWA uses 
another kind of sparks called Gaussian sparks, to further ensure 
diversity of the swarm. Compared to the explosion sparks, Gaus¬ 
sian sparks have global search ability. 

In the history of FWA, two typical Gaussian mutation op¬ 
erators were proposed. In FWA, the Gaussian mutation opera¬ 
tor focuses on the mutation ability for each firework alone by 
multiplying a random value under Gaussian distribution with the 
positions of the firework in the selected dimensions. Later on in 
EFWA, the introduced Gaussian mutation operator emphasizes the 
cooperative evolution in the fireworks swarm (see Fig.j2|. 


(0,0) 
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Fig. 2. Gaussian mutation operator in FWA and EFWA 
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2.3.1 Generating Gaussian sparks 

Alg.[3]gives some details about how Gaussian sparks are generated 
in EFWA. This algorithm is performed M g times, each time with 
a randomly selected firework Xi ( M g is a constant to control the 
number of Gaussian sparks, Gaussian(a, b ) denotes a random 
value from a normal distribution with expected value and variance 
set to a and b, respectively). 

Algorithm 3 - Generating one “Gaussian spark” in EFWA 
1: Initialize the location of the “Gaussian sparks”: Xi = Xi 
2: Set z k = round(rand(0,1)), k = 1,2, d 
3: Calculate offset displacement: e = Gaussian(0,1) 

4: for each dimension Xi, where z k == 1 do 
5: Xik = x ik + (A'bfc — Xik) * e, where Xf t is the position of 

the best firework found so far. 

6: if Xik is out of bounds then 

7: x ik = X^ in + rand(0,1) * (X^ ax - X^ in ) 

8 : end if 

9: end for 


2.3.2 Mapping out-of-bounds sparks 

When the location of a new Gaussian spark exceeds the search 
range in the fcth dimension, the mapping method is similar as the 
explosion sparks which exceed the search range. 

2.4 Selection Strategy 

To retain the information of the swarm and pass them to the next 
iteration, a subset of the whole population has to be selected. 
EFWA applies a computationally efficient selection method called 
Elitism-Random Selection (ERP, (45|). The optimum of the set 
( i.e. all fireworks, the explosion sparks, and the Gaussian sparks) 
will be selected firstly and the rest of individuals are selected 
randomly. Although this method is rather simple compared to 
e.g., distance based selection method in FWA (cf. (T9), (4§), our 
analysis in |24[ has shown that there is almost no difference in 
terms of convergence, final fitness and standard deviation between 
the distance based selection of conventional FWA and ERP. 

2.5 Comparison with Other Meta-heuristic Algorithms 

Since its introduction, FWA has shown great success for opti¬ 
mization problems, which indicates that some mechanisms of this 
algorithm are efficient. Here, we pick two famous algorithms, 
particle swarm optimization (PSO) and genetic algorithm (GA) 
for comparison. 

2.5.1 Reproduction Mechanism 

From the offspring generation view, we present the following 
definitions: “ cross-propagation ’ means that the offspring is gen¬ 
erated by at least two parents, while “ self-propagation” means the 
offspring is generated by only one parent. 

Thus, it can be seen that both cross-propagation and self¬ 
propagation are presented in GA, and there is no bijection from 
generation to generation. In PSO, the position of the son particle 
is based on the information of this particle’s history information 
and the global information of the best particle. Only the cross¬ 
propagation lies in PSO and there is a bijection from generation 
to generation. The reproduction mechanism in FWA is similar to 
GA. The explosion sparks are generated by the explosion operator 


performed by one firework and mutation sparks are generated 
by the mutation operator, i.e. both cross-propagation and self¬ 
propagation are presented in FWA and there is no bijection. 
However, different from the GA, the explosion operator in FWA 
plays the major role for optimization. 

2.5.2 Exploration and Exploitation Capacities 

Previous work in (47 1 has given the comparison between PSO 
and GA. For GA, the crossover operator is usually very effective. 
At the early search phase, the crossover operator can move the 
chromosome to a very far space in the solution space, which 
makes GA maintain high exploration capacity. While at the later 
phase, as the population of chromosomes converges, the crossover 
operator will not be able to move it to a far space as the two 
chromosomes usually have the similar structure, which makes the 
GA maintain high exploitation capacity. In addition, the mutation 
operator can move the chromosome to anywhere in the solution 
space which enhances the exploration capacity of GA. In PSO, 
each particle is attracted by current position, its previous best 
position and the global best position (of the entire population). 
This has similar function as the crossover operator in GA (47) , 
and finally the whole population of particles will converge to one 
position, which makes the swarm maintain a high exploitation 
capacity. The particles have to move one by one step from one 
place to another, which cannot jump with a quite high step to 
another place in one iteration (47). However, it seems that GA 
and FWA can do this by uding a big mutation rate and a big 
explosion amplitude, respectively. For FWA, the exploration and 
exploitation capacities are gained by taking different explosion 
amplitudes and explosion sparks number for each firework. A 
small explosion amplitude with large explosion sparks number 
can make the firework perform a comprehensive local search. To 
some extend, FWA is quite similar to GA, the smaller explosion 
amplitude has the similar effect as mutation operator with small 
mutation rate while the big amplitude has the similar effect as 
mutation operator with big mutation rate. 

2.5.3 Survival strategy 

In GA, the selection strategy is performed in each iteration to 
maintain and pass the information to the next iteration. The fittest 
candidate will be surely selected to the next iteration, while for the 
rest of candidates, they are selected with probability, which means 
that there is no bijection from generation to generation. For PSO, 
it does not utilize the selection operator, and there is bijection from 
generation to generation for each particle. In FWA, the selection 
operator is performed like GA in each iteration. Section |4.1| will 
give a detailed discussion with the survival strategy. 

3 State-of-the-art Improvement Work: 
dynFWA and AFWA 
3.1 Definition 

For the convenience of discussion, the following definitions are 
introduced at first: 

Core firework (CF) and non-core firework (non-CF) : Among the 
fireworks swarm, the firework with minimal fitness according to 
the / is denoted as CF. All other fireworks except the CF are 
denoted as non-CFs. 

General core firework (GCF) and non-general core firework 
(non-GCF) : The fireworks which have the same parent firework 
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as current CF may have the similar performance as the CF due to 
the selection method. Thus, we extend the CF to general CF by 
including the fireworks which have the same parent as CF. The 
detailed definition is given below. 

In each iteration, after the generation of explosion sparks, in 
the candidates set ( S c ), which includes the fireworks and explosion 
sparks, the selection operation is performed. Let us define x b as 
the “best” newly created explosion spark of all fireworks in the 
swarm, and Xcf as the current core firework. Assume that P{xi) 
denotes the firework which generated Xi, and Sqcfs denotes the 
set of GCF. 

• If f{x b ) — f(XcF ) < 0, and the candidate Xi which 
is from the candidate set S c is subjected to the following 
condition, 

Xi = arg (P(xi) == P(x b ))\\(xi == P(x b )) 

Xi 

and Xi is selected as firework into next iteration, then x^ £ 

Sgcfs- 

• If f{x b ) — J{X C f) > 0, and the candidate x, which 
is from the candidate set S c is subjected to the following 
condition, 

Xi = arg (P{xi) == X C F)\\(xi £ S G cf s ) 

Xi 

and Xi is selected as firework into next iteration, then Xi £ 

Sgcfs- 

For the non-GCF, among the fireworks swarm, the fireworks 
except for GCF are denoted as non-GCFs. 

3.2 Explosion Strategies in dynFWA and AFWA 

In EFWA, the MEACS (ref Section [272] ) enforces the exploration 
capabilities at the early phase of the algorithm (larger A m - ln => 
global search), while at the final phase of the algorithm the ex¬ 
ploitation capabilities are enforced (smaller A m ; n => local search). 
Additionally, the non-linear decrease of A m j n (cf. Eq. [4]) enhances 
exploitation already at an earlier stage of the EFWA algorithm. 
However, this procedure decreases the explosion amplitude solely 
with the current number of function evaluations which heavily 
depends on the pre-defined number of iterations for the algorithm. 
The explosion amplitude strategy should consider the optimization 
process information and the information of the explosion sparks 
rather than solely the information about the current iteration 
[evaluation] count. To approach this problem, the adaptive explo¬ 
sion amplitude strategies for fireworks were proposed, which are 
dynamic search in FWA (dynFWA) and adaptive FWA (AFWA). 

In dynFWA and AFWA, the fireworks are separated into two 
groups. The first group consists of the CF, while the second group 
consists of all remaining fireworks. The responsibility of the CF 
is to perform a local search around the best location found so far, 
while the responsibility of the second group is to maintain the 
global search ability. 

3.2.1 The dynFWA 

Let us define x b as the “best” newly created explosion spark of all 
fireworks in the swarm, for the dynFWA, the explosion amplitude 
adapts itself according to the quality of the generated sparks. 

Case 1) : One or several explosion sparks have found a better 
position. It is possible that ( i) an explosion spark generated by 
the CF has found the best position, or that (ii) an explosion spark 


Algorithm 4 - dynFWA explosion amplitude update strategy 
Require: Define: 

Xcf is the current location of the CF; 

x b is the best location among all explosion sparks; 

A * l CF is the CF’s explosion amplitude in iteration t; 

C a is the amplification factor; 

C r is the reduction factor; 

Ensure: 

l: if /(x b ) - f{X C F) < 0 then 
2: A^p <— A^p X C a \ 

3: else 

4: Agp <— J^Qp X C r \ 

5: end if 


generated by a different firework than the CF has found the best 
position. Both cases indicate that the swarm has found a new 
promising position and that x b will be the CF for the next iteration. 

(i) In most cases, x b has been created by the CF. In such 
cases, in order to speed up the convergence of the algorithm, 
the explosion amplitude of the CF for the next iteration will be 
increased compared to the current iteration. 

(ii) In other cases, a firework different from the CF will create 
x b . In such cases, x b will become the new CF for the next iteration. 
Since the position of the CF is changed, the current explosion 
amplitude of the current CF will not be effective to the newly 
selected CF. However, it is possible that x b is located in rather 
close proximity to the previous CF due to the fact that the random 
selection method may select several sparks created by the CF, 
which are initially located in close proximity to the CF. If so, the 
same consideration as in (i) applies. If x b is created by a firework 
which is not in close proximity to the CF, the explosion amplitude 
can be re-initialized to the pre-defined value. However, since it 
is difficult to define “close” proximity, we do not compute the 
distance between x b and Xcf but rely on the dynamic explosion 
amplitude update ability. Similarly to (i), the explosion amplitude 
is increased. If the new CF cannot improve its location in the next 
iteration, the new CF is able to adjust the explosion amplitude 
itself dynamically. 

We underline the idea why an increasing explosion amplitude 
may accelerate the convergence speed: Assume that the current 
position of the CF is far away from the global/local minimum. In¬ 
creasing the explosion amplitude is a direct and effective approach 
in order to increase the step-size towards the global/local optimum 
in each iteration, i.e., it allows for faster movements towards the 
optimum. However, we should make sure that the fireworks walk 
towards to global/local optimal position. In fact, the probability 
to find a position with better fitness decreases with increasing 
explosion amplitude due to the increased search space, it will be a 
hard job and a challenge for fireworks. 

Case 2) : None of the explosion sparks of the CF nor of all other 
fireworks has found a position with better fitness compared to 
the CF, i.e., A f >= 0. Under this circumstance, the explosion 
amplitude of the CF is reduced in order to narrow down the search 
to a smaller region around the current location and to enhance 
the exploitation capability of the CF. The probability for finding 
a position with better fitness usually increases with the decreasing 
of explosion amplitude. 
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3.2.2 The AFWA 

In AFWA, the motivation is to calculate the most suitable ampli¬ 
tude for the CF without preseted parameters by making full use of 
the region’s information based on the generated explosion sparks 
around the current position |26| . 

AFWA aims at finding a specific spark and using its distance 
to the best individual (which is either the Xcf or Xb , and will be 
selected as the CF for the next iteration) as the explosion amplitude 
for the next round. The specific spark that AFWA chooses should 
be subject to the following two conditions: 

• Its fitness is worse than the one of the current CF, which 
guarantees that the specific spark is not too close to CF; 

• The distance between the specific spark with the best 
individual is minimal among all the candidates, which is 
to ensure the convergence. 

As in FWA, the explosion sparks are generated independently 
in each dimension, the infinite norm is taken as the distance 
measure as follows. 


ll x lloo= maX (W) ( 5 ) 

l 

The AFWA’s explosion amplitude update strategy includes two 
cases: 

Case 1) The CF does not generate sparks with smaller fitness than 
the firework’s, i.e. the best fitness of all the explosion sparks XcF,b 
is worse than the core firework’s fitness Xcf, and the explosion 
amplitude will be set to the infinite norm between the core firework 
and the specific spark. If so, the explosion amplitude in the next 
iteration will be reduced according to Alg. [5] 

Case 2) The generated XcF,b has smaller fitness than Xqf, 
according to Alg. [5] it will be hard to determinate whether the 
explosion amplitude will be increased or decreased. Generally, the 
explosion amplitude will be increased with high probability while 
be reduced with smaller probability. 

Moreover, as the infinite norm between the specific spark and 
the firework changes radically, AFWA introduces the smoothing 
strategy (A^p 0.5 • (Agp + A • Aqp)). 

Algorithm 5 AFWA explosion amplitude update strategy 
Require: Define: 

M is explosion sparks number of CF; 

Aqp is the CF’s explosion amplitude in iteration t; 

XcF.k is the k- th explosion spark of CF; 

XcF.b is the best explosion spark of CF; 

A is a constant parameter which controls the update rate. 
Ensure: 

l: for k = 1 to M do 

2: if | \xcF,k —xcf, b I loo > Aqp tind f(xcF,k) > K^cf) 

then 

3: A^p «— || XcF,k — Xcf, fclloo 

4: end if 

5: end for 

6: Aqf 0.5 • (A^p + A • A^p) 


3.3 dynFWA vs AFWA 

From Alg. [4] and Alg. [5] it can be seen that if the CF cannot gen¬ 
erate any sparks with better fitness, then the explosion amplitude 
of the firework will be reduced in both AFWA and dynFWA. If 



(a) dynFWA, the explosion amplitude (b) dynFWA, the explosion amplitude 
is updated when no better solutions is updated when a better solution is 
are found. found. 




(c) AFWA, the explosion amplitude is (d) AFWA, the explosion amplitude 
updated when no better solutions are is updated when a better solution is 
found. found. 


Fig. 3. The comparison of explosion amplitude update between AFWA 
and dynFWA. 


the CF has generated a spark with better fitness, the explosion 
amplitude of the firework in dynFWA will be amplified. However, 
for AFWA, the determination of amplification or reduction is more 
complex. 

To investigate the relationship between of the amplification 
and reduction with whether a better solution is found or not, 
we record the values of Afp/Apfp. For dynFWA, if a better 


solution is found, the recorded value Af np /A^ C p will be C a 


^CF/^CF- 

^CF/^CF 

and the explosion amplitude is updated with the amplification 
factor. Otherwise, the recorded value A^p/A^p will be C r and 
the explosion amplitude is updated with the reduction factor (cf. 
Fig. |3(b)| and Fig. |3(a)| >. 

To investigate the explosion amplitude update process of 
AFWA, Sphere function is chosen as the benchmark function to 
study the explosion amplitude amplification and reduction when 
the fire work has found or not found a better solution. 

Fig. 3(d) presents the histogram of A^p/A^p in the iterations 
when a better s olution is found, here the y-coordinate denotes the 
probability. Fig. 3(c) presents the histogram of Af CF /Aeffp in the 
iterations when no better solutions are found. The geometric mean 
in Fig. |3(d)| is 1.017, while the geometric mean in Fig. |3(c)| is 
0.976. It can been seen that if the firework has found a better 
position with smaller fitness, the explosion amplitude in AFWA 
will be increased with higher probability, and reduced with smaller 
probability. 

Here, a special case is also considered to compare the two 
explosion amplitude strategies. Assume that explosion sparks 
number M is a very big value, i.e. M -A +oo, and the gen- 
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Fig. 4. dynFWA vs AFWA in a special case 


erated sparks are located in a float and continues region. Under 
this circumstance, if the generated sparks have found any better 
positions, usually the best explosion spark will be located at the 
edge of the region (see Fig. [4] with dimension is set to 2), then 
both AFWA and dynFWA will increase the explosion amplitude 
for the next iteration according to Alg[4]and Alg|5] 

Thus, we can draw the following conclusions: 

• The explosion amplitude is one of the most crucial pa¬ 
rameters for FWA to balance the global and local search 
ability. 

• From the statistical view, the two methods - dynFWA 
and AFWA are from different start points to adjust the 
explosion amplitude, but in a sense, finally they reach the 
same goal by different means. 


4 The Proposed FWA Cooperative Frame¬ 
work 

In this section, we first give a detailed analysis of the cooperative 
strategies in conventional FWA framework and the limitations 
are pointed out, and then the new cooperative FWA framework 
(CoFFWA) is finally proposed, which includes the independent 
selection method for each firework and crowdness-avoiding coop¬ 
erative strategy among the fireworks. 

4.1 Analysis of Cooperative Strategies in Conventional 
FWA Framework 

FWA is designed as one swarm intelligence algorithm, in which 
the fireworks in the swarm cooperate with each other to deal with 
the task that one firework cannot work well on. In EFWA, two 
cooperative operations are used for the implementation of this 
idea: 


• the sharing of the fitness values in the population for the 
calculation of explosion amplitudes and explosion sparks 
number in the explosion operator (cf subsection |2.2[ > ; 

• the Gaussian mutation strategy in the Gaussian mutation 
operator (cf subsection |2.3| >. 

The sharing of fitness values among the fireworks to calculate 
the explosion amplitudes and sparks numbers makes the fireworks 
with smaller fitness have smaller explosion amplitudes and larger 
number of explosion sparks, which maintains the capacity of 


exploitation while the fireworks with bigger fitness have bigger 
explosion amplitudes and smaller numbers of explosion sparks, 
which maintains the capacity of exploration. As for the Gaussian 
mutation operator, the generated Gaussian mutation sparks will 
locate along the direction between the core firework and the 
selected firework which is assumed to improve the diversity of 
the fireworks swarm. 

However, for EFWA, dynFWA, and AFWA, all of them take 
the selection method from the candidates set which includes the 
fireworks in current iteration, the generated explosion sparks and 
Gaussian mutation sparks in each iteration. For all of the FWA 
variants, the basic principle of the selection method is that the 
optimum among the candidates is always kept while for the rest 
of fireworks, different methods have different selection probability 
calculation methods. 

In the following, we will present an observation that these 
kinds of selection methods result in the fact that in the explosion 
operator the non-CFs contribute less for the optimization while 
taking up the most evaluation times, and the Gaussian mutation 
operator proposed in EFWA is not as effective as it was designed 
to be. 

4.1.1 Cooperative Performance Analysis of CF and non- 
CFs in Explosion Operator 

In FWA, the calculation method of explosion amplitude will make 
the explosion amplitude of the CF close to zero. To avoid this 
limitation, EFWA introduces the MEACS (cf Section |2.2[ >, where 
the CF’s explosion amplitude is in fact calculated based on the 
MEACS which is not relevant to the non-CFs in the population. 
In dynFWA or AFWA, the situation is similar, that the explosion 
amplitude of CF is calculated based on the dynamic search strategy 
or adaptive strategy which is also not relevant to the non-CF’s 
fitness. Compare the CF in EFWA, dynFWA and AFWA, it can be 
seen that the CF’s explosion amplitude strategies are independent 
of the non-CFs’ in the fireworks swarm, the only interaction 
between the CF and non-CFs except for the selection method is 
the calculation of explosion sparks number which is powerless for 
dealing with complex problems. 

Comparing the CF with non-CFs, the differences between 
them lie in two aspects, one is that the explosion amplitudes’ 
calculation methods, the other is the selection method. For the 
explosion amplitude calculation method, different explosion am¬ 
plitude calculation methods are taken to maintain the diversity of 
the fireworks swarm. However, the selection method makes the 
initial idea that the non-CFs are to increase the diversity of the 
fireworks swarm not effective. 

In the selection, most of the sparks generated by CF and non- 
CFs are brought together. The candidate with minimal fitness [will 
be the CF in the next iteration] is always selected at first, while 
for the rest of fireworks [will be non-CFs in the next iteration], 
they are selected with probability. If a selected non-CF candidate 
is generated by the CF, then it will have similar performance 
as the CF. Alternatively, the selected candidate is generated by 
the non-CFs. In this case, the non-CFs are usually generated by 
the parent fireworks within large explosion explosion amplitudes 
which can be seen as randomly generating of sparks within the 
search range as the information of the non-CF can only be passed 
to the next iteration within few iterations with the conventional 
selection method, thus after a number of iterations, the position of 
the non-CF will be reinitialized in a new position and will be kept 
only in few iterations. 
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In fact, these kind of selection strategies will make the gen¬ 
erated sparks by the non-CFs have the similar performance with 
randomly generated sparks in the search space to some extend or 
have the similar performance with sparks generated by CF. We 
cannot expect the good fitness by randomly generating the sparks 
if there is no much heuristic information. 

In Section [5] experiments are designed to evaluate the perfor¬ 
mance of CF and GCF to validate our idea that the non-CFs and 
non-GCFs make little contribution to the optimization. 

4.1.2 Cooperative Strategy in Gaussian Mutation Operator 

The motivation of Gaussian sparks is to increase the diversity 
of the swarm. The prerequisites for successful performance by 
cooperating in the fireworks population is that the information of 
each firework is different and effective. 

In EFWA, the newly generated sparks locate along the di¬ 
rection between the selected firework (Xfi and CF ( Xqf )■ For 
the selected firework Al,, it can be classified into two categories. 
The first category comprises fireworks which are very close to 
Xqf , usually these fireworks have the same parent firework as the 
CF. If so, the newly generated Gaussian sparks may have similar 
performance with the explosion sparks which cannot increase 
the diversity of the swarm effectively. The other category is the 
fireworks which are not close to the CF, usually it is because these 
fireworks are generated by a firework different from the parent 
firework of CF. If so, the newly generated Gaussian sparks will be 
(«) close to CF ( ii ) close to the selected firework, (Hi) not close 
to any fireworks. If the newly generated Gaussian spark is close to 
either of the CF or the selected firework, it has similar performance 
with the explosion sparks generated by them. If the generated 
Gaussian spark is not close to them, then it can be seen as the 
spark generated by firework with large explosion amplitude. Thus, 
based on above analysis, the newly generated Gaussian spark will 
not work effectively and is not be able to increase the diversity of 
the fireworks population. 

Moreover, assume that the selected firework X, is a non-CF 
with a random position due to the selection method, then the 
generated Gaussian mutation spark will have the similar effects 
with randomly generating a spark in the search range under the 
case (ii) and case (Hi). Thus, the generated Gaussian sparks will 
not be able to improve the diversity of the fireworks swarm. 

4.1.3 Conventional FWA Framework vs Evolutionary Strat¬ 
egy 

In the previous sections, we have pointed out that the non-CFs 
and Gaussian sparks do not contribute much to the optimization 
due to the selection method in the conventional FWA framework. 
If we simply eliminate the Gaussian mutation operator and non- 
CFs in the explosion operator, we will get the minimalist fireworks 
algorithm (MFWA) using only one firework, as shown in Alg. [6] 


Algorithm 6 The MFWA 


1: Initialize the firework and explosion amplitude 

2 : repeat 

3: Generate explosion sparks (cf. Alg. [2 i 

4: Update explosion amplitude (cf. Alg. 

5: Select the best individual as the new firework 

6: until termination criterion (time, max. # evals, convergence, 
... ) is met 


4] or Alg. [5} 


For comparison, we briefly describe a kind of evolutionary 
algorithms: evolutionary strategy (ES) |48| . Generally speaking, 
ES conducts the search process by keeping the iterations of 
sampling and learning. In the sampling step, a number of points 
are randomly sampled from a Gaussian distribution. In the learning 
step, the parameters of the Gaussian distribution are updated 
according to the qualities of the sampled points. 

The iteration process of MFWA is very similar to a (1 + A)- 
ES: one parent generates many sons, and the sampling parameters 
are updated according to the qualities of them (the position and 
the fitness). The explosion amplitude in FWA and the variance in 
ES can both be regarded as the step size. The methodology to 
control them are also very similar: dynamic explosion amplitude 
is similar to the 1/5 rule in (1 + 1)-ES, while adaptive explosion 
amplitude is comparable to (1 + A)-ES. So the properties of these 
algorithms are similar: they are all locally convergent, they all 
work poor in dimension sensitive functions, and they are all easily 
trapped in local minima. After the comparison with ES, it leaves a 
profound insight that the new FWA framework which can utilize 
all the fireworks information should be designed and the powerful 
and efficient cooperative mechanisms are essential for the future 
developments. 

4.2 Cooperative Framework of FWA 

To make FWA be a successful swarm intelligence algorithm [45) , 
the information for each firework should be passed to the next 
iteration and fireworks can cooperate with each other for the 
optimization. 

However, the conventional FWA framework lacks the local 
search ability for non-CFs, while the cooperative strategy in the 
Gaussian mutation operator is not very effective. To tackle these 
limitations, the CoFFWA with independent selection method and 
cooperative strategy is finally proposed. 


Algorithm 7 - The Cooperative Framework for FWA 
1: Initialize N fireworks 

2 : repeat 

3: for each firework do 

4: Generate the explosion sparks 

5: Select best candidate among its generated sparks and 

itself 

6: Perform the cooperative strategy 

7: end for 

8: until termination criterion (time, max. # evals, convergence, 
... ) is met 


4.2.1 Independent Selection Method 

The analyses in previous subsection suggest that the previous 
random selection methods in FWA and its variants cause that non- 
CFs and non-GCFs do not contribute much to the optimization 
for a problem while consuming a lot of resources, i.e. evaluation 
times in the explosion operator, and the cooperative scheme of 
Gaussian mutation operator among the fireworks is powerless to 
solve complex problems. Furthermore, the selection method in 
conventional FWA framework is seen as the main reason why 
the non-CFs/non-GCF cannot evolve for a number of consecutive 
iterations. 

To implement the initial idea of FWA that each firework in the 
swarm cooperate together to solve the optimization problem, it is 
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reinitilization 



Fig. 5. The Crowdness-avoiding Cooperative Strategy. 

needed to ensure that the information for each firework is passed to 
the next iteration. In the new framework of FWA, the independent 
selection operator will be performed for each firework, i.e. each 
firework will select the best candidate among its all generated 
sparks and itself in each iteration respectively (cf. Alg|7|. 

4.2.2 The Crowdness-avoiding Cooperative Strategy 
among the Fireworks 

In the fireworks swarm, each firework will generate a number 
of explosion sparks which can represent the quality of the local 
regions. The quality of each firework’s position is shared in the 
population to accelerate the convergence speed. 

For the calculation of explosion amplitudes and explosion 
sparks number in CoFFWA, CF still takes the dynamic explosion 
amplitude strategy, while for the rest of fireworks, the explosion 
amplitudes are calculated as in dynFWA by Eq. 0 - For explosion 
spark number calculation, it is same in CoFFWA as in dynFWA 
as well by Eq.Q. After the generating of explosion sparks, each 
firework performs the independent selection method, respectively. 

Moreover, In COFFWA, a crowdness-avoiding operation is in¬ 
troduced in the fireworks population. For the crowdness-avoiding 
operation, it means that whenever any fireworks in the fireworks 
swarm is close to CF, i.e. within a fixed range of CF, then the 
position of this firework will be reinitialized in the feasible search 
space. The following gives two main reasons for the crowdness- 
avoiding operation: 

• The selected firework ( Xi ) has worse fitness and smaller 
explosion sparks number than the CF, thus the region 
that Xi locates is not promising, it is seen as a waste of 
evaluation times to continue the search for the fireworks 
with worse fitness around the CF. 

• To increase the diversity of the fireworks swarm, it is better 
that the fireworks are not located close to each other. 


Algorithm 8 - The crowdness-avoiding cooperative strategy 

1: if | |.Xj - A'gfIIoo < r then 

2: Reinitialize the position of Xi 

3: end if 


5 Experiments Design 

To validate the ideas presented previously, several groups of 
experiments are designed. 


5.1 Significance Analysis of non-CFs/non-GCFs in Ex¬ 
plosion Operator 

To investigate whether the non-CFs/non-GCFs are effective or 
not compared with the CF/GCF in the explosion operator, the 
following evaluation criterions, i.e. significant improvement and 
resource cost are designed. 

For these following criterions evaluation, FWA variants with¬ 
out Gaussian mutation operator, i.e. the algorithms only generate 
explosion sparks are used to offset the influences of Gaussian 
mutation operators. 

5.1.1 Significant Improvement 

Among the fireworks swarm, if a firework Xi generates the 
explosion spark with minimal fitness among all the explosion 
sparks and fireworks, then Xi is seen to have one time significant 
improvement to the optimization. In the optimization process 
of FWA, at each iteration, at most one firework can make the 
significant improvement, thus to compare the performance of CF, 
GCF and non-CFs, non-GCFs, the significant improvement times 
of them during each run are recorded. 

• cxcf, the percentage of CF’s significant improvement 
times among all the significant improvement times. 

• /3c f> the percentage of CF’s significant improvement 
times among all the significant improvement times 
recorded from the A * .E max -th evaluation time. 

• &GCF, the percentage of GCF’s significant improvement 
times among all the significant improvement times. 

• Pgcf, the percentage of GCF’s significant improve¬ 
ment times among all the significant improvement times 
recorded from the ^ * i? max -th evaluation time. 

Here, E max denotes the max number of evaluation times. It 
is thought that /3cfI/3gcf are better than cxcfIoigcf for the 
performance comparison as at the early phase of the optimization, 
it is likely that non-CFs gain more significant improvement times 
due to the great number of explosion sparks. However, significant 
improvements in the later phase of the optimization are more 
important. 

5.1.2 Resources Cost 

For the optimization, the evaluation times (E max ) is usually set to 
D x 10 000, where D is the dimension of the problem )49| . 

• OcFi the percentage of CF’s evaluation times among all 
the evaluation times. 

• Oqcf• the percentage of GCF’s evaluation times among 
all the evaluation times. 

5.2 Significance Analysis of Gaussian Mutation Opera¬ 
tor 

To validate whether the Gaussian mutation operator introduced in 
EFWA is effective or not, the following comparisons are made. 

. EFWA-G vs EFWA-NG. 

• dynFWA-G vs dynFWA-NG. 

. AFWA-G vs AFWA-NG, 

here, “G” and “NG” refer to with and without Gaussian mutation 
operator, respectively. 
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5.3 Significance Analysis of the Proposed CoFFWA 

To validate the performance of the proposed CoFFWA, experi¬ 
ments are conducted to compare the performance of CoFFWA 
with EFWA (24), dynFWA (25) , AFWA |26|, artificial bee colony 
(ABC) (16), (50), differential evolution (DE) (5l|, standard PSO 
in 2007 (SPS02007) (4), and the most recent SPSO in 2011 
(SPSO2011) (52). 

5.4 Experimental Setup and Platform 

Similar to dynFWA, the number of fireworks in CoFFWA is set to 

5 and the explosion sparks number is set to 150. The reduction and 
amplification factors C r and C a are also empirically set to 0.9 and 
1.2, respectively. For the parameter r is set to 10*Acfi and the 
maximum explosion amplitude for CF is bounded with the search 
range. For the rest of parameters in CoFFWA, they are identical 
to dynFWA (25) . The parameters for AFWA, dynFWA, DE, ABC, 
SPS02007 and SPSO2011 are listed in (26) (16) (25) (5T) (4) (53), 
respectively. 

In the experiments, the CEC2013 benchmark functions with 
28 functions are used as the test suite and dimension is set to 
30 (49) . For the recorded results in each algorithm, the number 
of runs is set to 51, and Wilcoxon signed-rank test is used to 
validate the performance improvement. The experimental platform 
is MATLAB 2011b (Windows 7; Intel Core i7-2600 CPU @ 3.7 
GHZ; 8 GB RAM). 

6 Experimental Results 

6.1 Significance Analysis of non-CFs/non-GCFs in Ex¬ 
plosion Operator 

To validate the performance of CF and GCF, the acF- Pcf, 
9c f and accF > P GCF, 9qcf are calculated on EFWA-NG 
and dynFWA-NG respectively, and the recorded results on 28 
functions can be found in Fig. [6] 

Compare the performance between CF and non-CFs, it can be 
seen that, for both EFWA-NG and dynFWA-NG, CF takes smaller 
percentage of resources ( i.e. evaluation times), but makes the more 
significant improvement times to the search. 

Compare the evaluation criterion acF and Pcf, it can be 
seen that for all functions, Pcf which records the statistical results 
from the * .E max -th evaluation times is higher than acF, which 
means that at the beginning of the optimization, the non-CFs have 
high probability for making significant improvements, while at the 
later phase, the chance goes smaller. For GCF, the above situation 
is similar to CF. 

Compare the results of CF with GCF, the GCF makes more 
significant improvements while taking more resources (i.e evalua¬ 
tion times), due to that some fireworks except for CF may be close 
to CF, thus to have a high chance to make improvements. 

In all. from the results in Fig. [6] it can be concluded that the CF 
makes much more contributions than non-CF while taking smaller 
resources. While comparing the results of CF and GCF, it seems 
that the non-GCFs do not make any contribution to the search. 

6.2 Significance Analysis of Gaussian Mutation Opera¬ 
tor 

Tab. |T| gives the experimental results of the versions of EFWA, 
dynFWA, AFWA with and without Gaussian mutation operator. 
Comparing EFWA with EFWA-NG, it can be seen that EFWA-G 


is significantly better than EFWA-NG on 4 functions, which sug¬ 
gests that Gaussian mutation operator is effective to improve the 
performance for EFWA. However, for dynFWA-NG and AFWA- 
NG, the versions without Gaussian mutation operator have better 
performance. Moreover, in terms of the computational complexity, 
our previous results in (25) suggest that Gaussian mutation opera¬ 
tor is more time consuming than explosion sparks operator. Thus, 
the Gaussian mutation operator will be removed from CoFFWA. 

TABLE 1 

Wilcoxon signed-rank test results for EFWA-G vs EFWA-NG and 
dynFWA-G vs dynFWA-NG and AFWA-G vs AFWA-NG (bold values 
indicate the performance difference is significant, while 1/0/-1 denotes 
the version with Gaussian sparks operator is significant better / not 
significant different / significant worse than the version without 
Gaussian sparks.). 


F. 

EFWA-G vs 
EFWA-NG 

dynFWA-G vs 
dynFWA-NG 

AFWA-G vs 
AFWA-NG 

1 

2.316E-03 

1 

1.000E+00 

0 

1.000E+00 

0 

2 

4.256E-01 

0 

9.328E-01 

0 

4.647E-01 

0 

3 

8.956E-01 

0 

2.339E-01 

0 

7.191E-02 

0 

4 

7.858E-01 

0 

7.492E-02 

0 

5.689E-03 

-1 

5 

4.290E-02 

1 

7.646E-02 

0 

5.239E-01 

0 

6 

1.654E-01 

0 

7.858E-01 

0 

9.030E-01 

0 

7 

9.552E-01 

0 

6.869E-01 

0 

2.728E-01 

0 

8 

9.776E-01 

0 

4.704E-01 

0 

8.808E-01 

0 

9 

5.178E-01 

0 

4.997E-01 

0 

7.571E-01 

0 

10 

3.732E-01 

0 

5.057E-01 

0 

9.545E-03 

-1 

11 

5.830E-02 

0 

6.629E-01 

0 

3.204E-01 

0 

12 

6.193E-01 

0 

3.783E-01 

0 

1.801E-01 

0 

13 

8.220E-01 

0 

1.863E-01 

0 

4.590E-01 

0 

14 

4.101E-02 

0 

3.834E-01 

0 

5.239E-01 

0 

15 

6.869E-01 

0 

3.438E-01 

0 

4.879E-01 

0 

16 

2.81 IE-01 

0 

4.256E-01 

0 

8.220E-01 

0 

17 

9.179E-01 

0 

1.863E-01 

0 

2.339E-01 

0 

18 

6.938E-01 

0 

4.762E-01 

0 

6.460E-01 

0 

19 

9.402E-01 

0 

1.542E-01 

0 

7.786E-01 

0 

20 

1.559E-02 

-1 

5.830E-02 

0 

4.997E-01 

0 

21 

6.910E-04 

1 

7.997E-01 

0 

5.937E-01 

0 

22 

9.776E-01 

0 

3.583E-01 

0 

4.202E-01 

0 

23 

7.217E-01 

0 

7.642E-01 

0 

6.260E-01 

0 

24 

1.079E-02 

1 

5.486E-01 

0 

7.500E-01 

0 

25 

8.734E-01 

0 

2.091E-01 

0 

1.369E-02 

-1 

26 

2.687E-01 

0 

7.217E-01 

0 

3.834E-01 

0 

27 

3.534E-01 

0 

8.734E-01 

0 

1.597E-01 

0 

28 

6.460E-01 

0 

0.000E+00 

-1 

3.285E-03 

-1 


6.3 Significance Analysis of CoFFWA 

Tab. [2] shows the mean fitness value over 51 runs for the 28 
functions for ABC, DE, SPS02007, SPSO2011, EFWA, AFWA, 
dynFWA and CoFFWA and the corresponding rank of each al¬ 
gorithm. It can be seen that CoFFWA is able to achieve better 
results than dynFWA(AFWA) on 24(22) functions while dyn- 
FWA(AFWA) is better than CoFFWA on 3(5) functions. For 1(1) 
function, the results are identical. For the comparison with EFWA, 
the advantage is especially obvious. 

The Wilcoxon signed-rank test results in Tab. [3] suggest that 
CoFFWA is significantly better than dynFWA(AFWA) on 8(8) 
functions while significant worse on 2(2) functions. Furthermore, 
in terms of computation complexity (cf Tab. |4), CoFFWA has the 
similar computation cost as dynFWA. Runtime results in Tab. [4] 
suggest that compared to EFWA, dynFWA, the proposed CoFFWA 
has almost the same runtime. 

Compared with ABC, DE, SPS02007, SPSO2011, EFWA, 
AFWA and dynFWA in term of mean fitness rank, it can be 
seen that CoFFWA achieves the best optimization results and 
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(b) Case 2: EFWA-NG, GCF 



function id 



(c) Case 3: dynFWA-NG, CF 


(d) Case 4: dynFWA-NG, GCF 


Fig. 6. Significant improvement and resources cost results of CF and GCF. 


performances, i.e. the minimal average of mean fitness rank (3.00 
in Tab. |2j, thus we can conclude that the proposed cooperative 
framework for FWA is significant. 


TABLE 3 

Wilcoxon signed-rank test results for CoFFWA vs AFWA and CoFFWA 
vs dynFWA (1/0/-1 denotes the CoFFWA is significant better / not 
significant different / significant worse than EFWA/AFWA/dynFWA). 


F. 

CoFFWA vs 

AFWA 

CoFFWA vs 
dynFWA 

1 

1.000E+00 

0 

1.000E+00 

0 

2 

6.938E-01 

0 

6.460E-01 

0 

3 

1.313E-01 

0 

9.328E-01 

0 

4 

0.000E+00 

-1 

0.000E+00 

-1 

5 

0.000E+00 

-1 

0.000E+00 

-1 

6 

2.160E-01 

0 

3.989E-01 

0 

7 

1.132E-01 

0 

3.889E-03 

1 

8 

5.802E-01 

0 

7.077E-01 

0 

9 

4.202E-01 

0 

9.850E-01 

0 

10 

2.332E-01 

0 

2.981E-01 

0 

11 

2.853E-01 

0 

2.413E-01 

0 

12 

2.528E-01 

0 

7.963E-02 

0 

13 

4.041E-01 

0 

5.867E-01 

0 

14 

1.334E-02 

1 

8.676E-03 

1 

15 

8.600E-05 

1 

8.900E-05 

1 

16 

3.486E-01 

0 

1.153E-01 

0 

17 

1.000E-06 

1 

3.000E-06 

1 

18 

2.091E-01 

0 

4.148E-01 

0 

19 

4.094E-01 

0 

2.567E-01 

0 

20 

1.132E-01 

0 

6.260E-01 

0 

21 

0.000E+00 

1 

1.000E-06 

1 

22 

2.489E-01 

0 

8.513E-01 

0 

23 

2.811E-01 

0 

1.197E-01 

0 

24 

8.440E-01 

0 

2.528E-01 

0 

25 

3.659E-02 

1 

3.783E-01 

0 

26 

3.300E-05 

1 

2.316E-03 

1 

27 

2.966E-02 

1 

5.528E-03 

1 

28 

2.865E-02 

1 

1.594E-03 

1 


8+/18/2- 


8+/18/2- 



7 Conclusion 

In this paper, we have presented a cooperative framework for FWA 
(CoFFWA). The contributions include three aspects. 

1) The cooperative strategies in conventional FWA are analyzed 
and evaluated. In the explosion operator, the performance 


TABLE 4 

Runtime Comparison 

EFWA | 1.3 || AFWA | — || dynFWA | 1 || CoFFWA | T05 


of CF, GCF, non-CFs and non-GCFs in EFWA, dynFWA 
and AFWA are investigated and the designed criterions are 
recorded to support our opinion that the non-CF/non-GCFs 
do not make much contribution to the optimization for a 
specific problem while taking a lot of evaluation times. For 
the Gaussian mutation operator, experimental results suggest 
that it is not as effective as it is designed to be. 

2) Based on the analysis of cooperative strategies in con¬ 
ventional FWA framework, a cooperative framework of 
FWA with an independent selection method and crowdness- 
avoiding cooperative strategy is finally proposed to ensure 
the information inheritance and improve the diversity of 
fireworks population. 

3) Moreover, the explosion amplitude strategies in AFWA and 
dynFWA are compared and the relationship between them is 
presented. 

Compared with PSO’s search manner— one individual gen¬ 
erates one individual , CoFFWA presents a cooperative explosive 
search manner — one individual generates a number of individu¬ 
als which has the ability to estimate the property of typical local 
search space and will make the search direction of the swarm 
stable and the cooperative strategies among the fireworks make the 
firework with better fitness will have more resources. In the future, 
we need to focus on this unique search manner and try to develop 
new kinds of cooperative search strategies among the fireworks 
swarm to enhance the exploration and exploitation capacities. 
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TABLE 2 

Mean Fitness and mean fitness rank of ABC, DE, SPS02007, SPSO2011, EFWA, AFWA, dynFWA and CoFFWA (AR denotes the average of 

mean rank value). 


F 

ABC 


DE 


SPS02007 


SPSO2011 


EFWA 


AFWA 


dynFWA 


CoFFWA 


1 

0.00E+00 

1 

1.89E-03 

7 

0.00E+00 

1 

0.00E+00 

1 

8.50E-02 

8 

0.00E+00 

1 

O.OOE+OO 

1 

O.OOE+OO 

1 

2 

6.20E+06 

8 

5.52E+04 

1 

6.08E+06 

7 

3.38E+05 

2 

5.85E+05 

3 

8.92E+05 

6 

8.71E+05 

4 

8.80E+05 

5 

3 

5.74E+08 

7 

2.16E+06 

1 

6.63E+08 

8 

2.88E+08 

6 

1.16E+08 

3 

1.26E+08 

5 

1.23E+08 

4 

8.04E+07 

2 

4 

8.75E+04 

7 

1.32E-01 

1 

1.03E+05 

8 

3.86E+04 

6 

1.22E+00 

2 

1.14E+01 

4 

1.04E+01 

3 

2.01E+03 

5 

5 

0.00E+00 

1 

2.48E-03 

7 

O.OOE+OO 

2 

5.42E-04 

3 

8.05E-02 

8 

6.00E-04 

5 

5.51E-04 

4 

7.41E-04 

6 

6 

1.46E+01 

2 

7.82E+00 

1 

2.52E+01 

4 

3.79E+01 

8 

3.22E+01 

7 

2.99E+01 

5 

3.01E+01 

6 

2.47E+01 

3 

7 

1.25E+02 

7 

4.89E+01 

1 

1.13E+02 

6 

8.79E+01 

2 

1.44E+02 

8 

9.19E+01 

4 

9.99E+01 

5 

8.99E+01 

3 

8 

2.09E+01 

6 

2.09E+01 

2 

2.10E+01 

7 

2.09E+01 

5 

2.10E+01 

8 

2.09E+01 

4 

2.09E+01 

3 

2.09E+01 

1 

9 

3.01E+01 

8 

1.59E+01 

1 

2.93E+01 

6 

2.88E+01 

5 

2.98E+01 

7 

2.48E+01 

4 

2.41E+01 

3 

2.40E+01 

2 

10 

2.27E-01 

5 

3.24E-02 

1 

2.38E-01 

6 

3.40E-01 

7 

8.48E-01 

8 

4.73E-02 

3 

4.81E-02 

4 

4.10E-02 

2 

11 

0.00E+00 

1 

7.88E+01 

3 

6.26E+01 

2 

1.05E+02 

7 

2.79E+02 

8 

1.05E+02 

6 

1.04E+02 

5 

9.90E+01 

4 

12 

3.19E+02 

7 

8.14E+01 

1 

1.15E+02 

3 

1.04E+02 

2 

4.06E+02 

8 

1.52E+02 

5 

1.58E+02 

6 

1.40E+02 

4 

13 

3.29E+02 

7 

1.61E+02 

1 

1.79E+02 

2 

1.94E+02 

3 

3.51E+02 

8 

2.36E+02 

4 

2.54E+02 

6 

2.50E+02 

5 

14 

3.58E-01 

1 

2.38E+03 

3 

1.59E+03 

2 

3.99E+03 

7 

4.02E+03 

8 

2.97E+03 

5 

3.02E+03 

6 

2.70E+03 

4 

15 

3.88E+03 

4 

5.19E+03 

8 

4.31E+03 

7 

3.81E+03 

3 

4.28E+03 

6 

3.81E+03 

2 

3.92E+03 

5 

3.37E+03 

1 

16 

1.07E+00 

5 

1.97E+00 

8 

1.27E+00 

6 

1.31E+00 

7 

5.75E-01 

3 

4.97E-01 

2 

5.80E-01 

4 

4.56E-01 

1 

17 

3.04E+01 

1 

9.29E+01 

2 

9.98E+01 

3 

1.16E+02 

5 

2.17E+02 

8 

1.45E+02 

7 

1.43E+02 

6 

1.10E+02 

4 

18 

3.04E+02 

8 

2.34E+02 

7 

1.80E+02 

4 

1.21E+02 

1 

1.72E+02 

2 

1.75E+02 

3 

1.88E+02 

6 

1.80E+02 

5 

19 

2.62E-01 

1 

4.51E+00 

2 

6.48E+00 

3 

9.51E+00 

7 

1.24E+01 

8 

6.92E+00 

5 

7.26E+00 

6 

6.51E+00 

4 

20 

1.44E+01 

6 

1.43E+01 

5 

1.50E+01 

8 

1.35E+01 

4 

1.45E+01 

7 

1.30E+01 

1 

1.33E+01 

3 

1.32E+01 

2 

21 

1.65E+02 

1 

3.20E+02 

6 

3.35E+02 

8 

3.09E+02 

3 

3.28E+02 

7 

3.16E+02 

5 

3.10E+02 

4 

2.06E+02 

2 

22 

2.41E+01 

1 

1.72E+03 

2 

2.98E+03 

3 

4.30E+03 

7 

5.15E+03 

8 

3.45E+03 

6 

3.33E+03 

5 

3.32E+03 

4 

23 

4.95E+03 

5 

5.28E+03 

6 

6.97E+03 

8 

4.83E+03 

4 

5.73E+03 

7 

4.70E+03 

2 

4.75E+03 

3 

4.47E+03 

1 

24 

2.90E+02 

7 

2.47E+02 

1 

2.90E+02 

6 

2.67E+02 

2 

3.05E+02 

8 

2.70E+02 

4 

2.73E+02 

5 

2.68E+02 

3 

25 

3.06E+02 

6 

2.80E+02 

1 

3.10E+02 

7 

2.99E+02 

5 

3.38E+02 

8 

2.99E+02 

4 

2.97E+02 

3 

2.94E+02 

2 

26 

2.01E+02 

1 

2.52E+02 

3 

2.57E+02 

4 

2.86E+02 

7 

3.02E+02 

8 

2.73E+02 

6 

2.61E+02 

5 

2.13E+02 

2 

27 

4.16E+02 

1 

7.64E+02 

2 

8.16E+02 

3 

1.00E+03 

7 

1.22E+03 

8 

9.72E+02 

5 

9.80E+02 

6 

8.71E+02 

4 

28 

2.58E+02 
AR: 4.14 

1 

4.02E+02 
AR: 3.18 

5 

6.92E+02 
AR: 5.04 

7 

4.01E+02 
AR: 4.64 

4 

1.23E+03 
AR: 6.79 

8 

4.37E+02 
AR: 4.25 

6 

2.96E+02 
AR: 4.43 

3 

2.84E+02 

AR: 3.00 

2 
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